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ABSTRACT 



in 



Context, HD 46375 is known to host a Saturn-like exoplanet orbiting at 0.04 AU from its host star. Stellar light reflected by the planet 
was tentatively identified in the 34-day CoRoT run acquired in October-November 2008. 

Aims. We constrain the properties of the magnetic field of HD 46375 based on spectropolarimetric observations with the NARVAL 
spectrograph at the Pic du Midi observatory. In addition, we use a high-resolution NARVAL flux spectrum to contrain the atmospheric 
parameters. With these constraints, we perform an asteroseismic analysis and modelling of HD 46375 using the frequencies extracted 
\ from the CoRoT light curve. 

Methods. We used Zeeman Doppler imaging to reconstruct the magnetic map of the stellar surface. In the spectroscopic analysis we 
fitted isolated lines using ID LTE atmosphere models. This analysis was used to constrain the effective temperature, surface gravity, 
and chemical composition of the star. To extract information about the p-mode oscillations, we used a technique based on the envelope 
autocorrelation function (EACF). 

Results. From the Zeeman Doppler imaging observations, we observe a magnetic field of ~ 5 gauss. From the spectral analysis, 
HD 46375 is inferred to be an unevolved K0 type star with high metallicity [Fe/H] = +0.39. Owing to the relative faintness of the 
star (rahip = 8.05), the signal-to-noise ratio is too low to identify individual modes. However, we measure the p-mode excess power 
and large separation Av = 153.0 ± 0.7 //Hz. 

Conclusions. We are able do constrain the fundamental parameters of the star thanks to spectrometric and seismic analyses. We 
conclude that HD 46375 is similar to a young version of a Cen B. This work is of special interest because of its combination of 
exoplanetary science and asteroseismology, which are the subjects of the current Kepler mission and the proposed Plato mission. 

Key words. Stars: planetary systems, Stars: oscillations, Stars: individual: HD 46375, Methods: data analysis, Methods: observational 
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,— ' , 1 . Introduction that the inclin ation of the planetary orbital plane is less than 83° 

■ (lHenry|2000l) . However, we have discovered that the light curve, 

t-H . The CNES CoRoT satellite has been orbiting the Earth at an al- which is the sum of me star emission and the light reflected by 

j> ; titude of 896 km, in a polar orbit, since late December 2006. the planeti is modu lated by the reflected planetary flux (Gaulme 

k> 1 ItS sclentlfic 8° als are dlvlded lnto two themes: asteroseismol- et al 2 010b, Paper II). This paper is dedicated to the analysis 

^ ■ ogy and the detection of transiting exoplanets. Apart from the of the stellar p rope rties, by means of the seismic analysis of the 

H ; 5-month long observation runs dedicated to the main targets, CoRoX data and spectro-polarimetric ground-based support. An 

. . , CoRoT performs short runs intended to observe a large set of analysis of the CoRoT light curve that is focused on the planet 

stars with different properties. We present the observation of w ju ^ e presented in Paper II 
the non-transiting extrasolar system HD 46375 in the astero- 

seismological field of CoRoT, during a 34-day run between HD 46375 is a solar-like star with properties similar to 
October 9 and November 1 1 2009. The objective is to obtain a Cen B that is particularly interesting because it may be a 
observational constraints on both the star and the planet, by per- younger version of the a Cen A+B system. It is an unevolved 
forming a seismic analysis of the star and detecting the changing K0 type star of high metallicity, with apparent magnitude from 
phases of the planet during its orbital period. We already know Hipparcos database m bip - 8.05. More over, from the rotation- 
that this planet does not transit in front of the star, which implies activity relationship dNoves et a l.l l 19841) and given the measure- 
ments of Ca ii H and K emission log R' H K = -4.94 and (B — V) 
' The CoRoT space mission, launched on 2006 December 27, was 86 reported in|Marcy et al.| (|2000|), we estimate the stellar rota- 
developed and is operated by the CNES, with participation of the tion period at 42+; days. The star is understood to be very quiet 
Science Programs of ESA, ESA's RSSD, Austria, Belgium, Brazil, with relative variations of the order of a few 100 ppm on time 
Germany and Spain. scales of a few days. 
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Fig. 1. Radial velocity associated with each exposure of the Sept- 
Oct 2008 NARVAL observations of HD 46375. The measure- 
ments are fol ded over the planeta ry orbital phase (3.02357 days), 
estimated bv lButler etal] d2006l) . 



Since early 2007, CoRoT has been observing solar-like oscil- 
lations in several solar-like stars, subgiants, and giants. However, 
the number of late-type main-sequence targets i s still small. 



Most of the main t argets have been F t ype sta rs (Barban et al 



2009; iGarcfa et al.1 l2009t iMosser et al.l l2009bt iBenomar et al 
2009). Photometric observations with CoRoT favour F-type 
stars, since they are typically brighter and also expected to ex- 
hibit higher amplitude oscillations. Asteroseismic spectroscopic 
ground-based campaigns are limited to bright stars (V < 6), and 
some G and early K type stars h ave been observed, e .g. i Hor 
(IVauclair et al.l 120081) . a Cen B dKieldsen et al.ll2005l) . t Ceti 
dTeixeira et akfl2009l) . and HD 49385 (Deheuvels et al. 2010, 



A&A, submitted). With HD 46375, we have the opportunity to 
investigate solar-like oscillations in a low-mass KO-type star, to 
improve the determination of its internal structure and to more 
tightly constrain the fundamental parameters of both the stellar 
and planetary system. 

To complement the CoRoT data, ground-based observations 
were carried out at the Pic du Midi observatory. We used the 
NARVAL spectropolarimeter to obtain both a high resolution 
visible spectrum and magnetic maps of the stellar surface, by 
means of Zeeman Doppler imaging. These observations allow 
us to monitor the changes in magnetic activity. 

We first present the analysis of the spectropolarimetric 
and spectrometric data (Sect. 2). We then present the CoRoT 
lightcurve (Sect. 3) and the seismic analysis (Sect. 4). We then 
combine the constraints to develop an asteroseismic model of 
HD 46375 (Sect. 5). Finally, we discuss the implications for the 
properties of the star-planet system (Sect. 6). 



2. Magnetic imaging and stellar parameters 

2.1. Spectropolarimetric observations 

The spectropolarimetric data were acquired with the NARVAL 
spectropolarimeter, permanently mounted at the 2-m Bernard 
Lyot Telescope at the Pic du Midi Observatory (France). 
Eighteen spectra were collected during a 35-night run between 
September 21 and October 27, 2008. A single spectrum cov- 
ers the full optical domain from 3700 to 10480 A at a spec- 
tral resolution of 65 000. The circular polarization is measured 
in a four-exposure sequence with different configurations of the 
wave plates. The data were pro cessed with the Libre-Esprit soft- 
ware (lDonatietal.ll 19971 12007h . 
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Fig. 2. Maximum-entropy fits (thin red line) to the Stokes V LSD 
profiles observed by TBL/NARVAL (thick black line) of HD 
46375. The rotational phase (in radians) of each observation and 
lcr error bars are also shown next to each profile. 

Table 1. Journal of observations. From left to right: heliocen- 
tric Julian date at mid-exposure (Col. 1); total exposure time 
(Col. 2); peak signal-to-noise ratio (per 2.6 km s _1 velocity bin) 
of each four-subexposure polarisation sequence (Col. 3); rms 
noise level (relative to the unpolarized continuum level I c and 
per 1.8 km s _1 velocity bin) in the circular polarization profile 
produced by least squares deconvolution (LSD) (Col. 4); orbital 
cycle (Col. 5); rotational cycle (Col. 6); and the radial velocity 
associated with each observation (Col. 7). 
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We calculated the stellar rotational and planetary orbital 
phase corresponding to our observations, by assuming a stellar 
rotation period of 42 days, as obtained from the rotati on-activity 
relationship, and the planetary orbital period from iButler et aH 
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Fig. 3. Maximum-entropy reconstruction of the large-scale mag- 
netic topology of HD 46375 as derived by the TBL/NARVAL 
data set. The radial, azimutal, and meridional components of the 
field are shown from top to bottom. Magnetic flux values are 
given in gauss as indicated along the vertical colour bar. The 
star is shown in flattened polar projection down to latitudes of 
-30°, with the equator depicted as a bold circle and parallels as 
dashed circles. Radial ticks around each plot indicate the rota- 
tional phases of observations. 

(2006 fl Least squares deconvolution (LSD) techniques were ap- 
plied to the spectra to imp rove the signal-to -noise ratio (SNR), 
with a K0 template mask dDonatiet al.lll997h . More than 6000 
stellar lines were used to derive mean profiles in the Stokes I and 
V parameters. The mean I profile was used to deri ve the radial 
veloci ty (Fig.|TJ, whose variation is compatible with lButler et alJ 
d2006l) . The error in the radial velocity (about 30 m s _1 per indi- 
vidual measurement) and the poor orbital sampling (with a pe- 
riod close to an integer number of days) do not allow us to refine 
the planetary ephemeris, since the semi-amplitude of the radial- 



HJD = 2454480.0 + 42.0 <p mt (1) 
HJD = 2451071.359 + 3.023573 <p mb (2) 



velocity variation is only 34 m s 1 . We present the journal of 
observations in Table [3] 

2.2. Magnetic imaging 

We use Zeeman Doppler imaging (ZDI) to reconstruct the mag- 
netic map of the stellar surface. ZDI is a tomographic imaging 
technique based on the inversion of the time se ries of the Stokes 
V profiles into a map of t he magnetic topology (Dona ti & Brownl 
ll997UDonati et al1l2006l) . An example of this modelling is given 
bv lFares et al.l d2009l) . who applied it to the planet-bearing star r 
Boo. 

The LSD Stokes V signature is detected in two thirds of the 
spectra, with small variations in the polarized profile (Fig. 0. 
Those variations over 35 nights are limited; in particular, the 
field component along the line of sight is always positive and 
does not switch sign across the rotation cycle, suggesting that 
the large-scale field topology is fairly simple and mostly axisym- 
metric with respect to the rotation axis of the star. We note that 
we checked the V Stokes profiles are significantly detected. Each 
profile shown in Fig. [2] has indeed a flat counterpart in the null 
profile, used for check, and no flat noise pattern was observed in 
the present data. 

The fact that the magnetic Stokes V measurements do not ex- 
hibit any periodicity, confirms that the rotation period is longer 
than 35 days. Therefore, the temporal coverage of the ZDI ob- 
servations is shorter than the expected rotation period. Hence, it 
is not possible to constrain either the rotation period or the dif- 
ferential rotation at the stellar surface. We use the 42^ day es- 
timate of the rotation period to invert the data. In addition, from 
the stellar modelling performed thanks to the seismic and spec- 
trometric analysis CSect. l5.3l l. the stellar radius is 0.91 + 0.02 R0. 
Hence, the rotational velocity at the equator is v = 1.1 ± 0.2 km 
s _1 . In Sect. 12. 31 the projected rotation velocity is estimated to be 
vsirn = 1.2 + 0.4 km s . This estimate is not fully compatible 
with the equatorial velocity estimate, whose upper limit is 1.3 
km s _1 this also being the v sin i upper limit. Hence, v sin i lies in 
the range [0.8, 1 .3] km s _1 . We then have a rough estimate of the 
stellar inclination: by assuming Gaussian distributions for v sin i 
and v, we obtain an estimate of the inclination angle i = 50+18°. 

For simplicity, we modelled the Stokes V profiles by taking 
an inclination of i = 45° (well within the error bars) and a line- 
of-sight projected equatorial rotation velocity v sin i = 0.9 km 
s _1 . The magnetic topology derived with ZDI is shown in Fig. [3] 
with the corresponding maximum-entropy fit to the data being 
shown in Fig. [2] The reconstructed field structure resembles a 
simple dipole slightly tilted, by about 15°, with respect to the ro- 
tation axis. The magnetic strengh at the pole is about 5 gauss. We 
note that the result is only moderately sensitive to the true value 
of i and changes by about +10° produce very small changes in 
the corresponding magnetic image. The stellar inclination is an 
input of the magnetic modelling and not a fitted output parame- 
ter. In other words, the magnetic mapping does not allow us to 
constrain the stellar inclination. 

The magnetic structure of HD 46375 appears to be signifi- 
cantly different from those of other planet-bearing stars such as 
t Boo dCatala et al J 120071: iDonati et alj|2008t iFares et al.l F2009) 
and HD 189733 dMoutou et al.ll2007l Fares et al., submitted to 
MNRAS). HD 46375 presents a far more simple magnetic struc- 
ture. More specifically, the magnetic field of these 2 stars shows 
a significant or even mainly toroidal component whereas that of 
HD 46375 is predominantly poloidal. This difference mostly re- 
flects the differences in rotation periods and more precisely the 
difference in Rossby numbers, which is related to the dynamo 
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prope rties for stars of different spectral types (e.g. Nov es et"a"D 
11984 . With a rotation period of 42 days, HD 46375 features a 
Rossby number of about 2 and naturally falls among stars featur- 
ing very simple, mo stly-poloidal and axisymme tric large-scale 
magnetic topologies dDonati & Land street 2009). 

2.3. Stellar parameters 

The high-resolution NARVAL spectra of HD 46375 were also 
used to determine the stellar parameters. All intensity spectra 
were summed prior to the spectroscopic analysis: the resulting 
SNR is 330 around 6000 A. The summed spectrum was anal- 
ysed in the spectral range 5000 to 8500 A with the VWA soft- 
ware, which make s use of iterative fitting of synthetic profiles 
dBruntt et alj|2004l 120081) . In Fig. |4] a portion of the observed 
spectrum is compared with its corresponding synthetic profile 
(dash-dotted line). We note that the star is late type (K0) and that 
several weak lines are not available i n the adopted atomic line 
data from VALD (Kup ka et al.l [l999). This introduces a weak 
bias in the abundance estimates. It may also affect the estimate 
of the atmospheric parameters, since weaker lines are more con- 
taminated by unidentified lines than relatively stronger lines. 

We used 159 isolated Fe i lines to constrain the atmospheric 
parameters, by minimizing the correlations between abundance 
with equivalent width and excitation potential. The neutral and 
ionized Fe lines were required to infer the same abundance, to 
constrain the surface gravity to be logg = 4.73 ± 0.10 dex. 
Furthermore, to confirm the \ogg value, the synthetic profiles 
were compared to the Mgib lines and the broad Ca lines at 
6122, 6262, and 6439 A dBruntt et al.ll2010l) . On the one hand, 
the Mg ib lines could not be used because this region is strongly 
blended and affected by molecular lines for such a late-type star. 
On the other hand, the determined surface gravity with the three 
Ca lines, is, respectively, logg = 4.64 ± 0.07 dex, 4.63 ± 0.04 
dex, and 4.69 ± 0.16 dex. These results agree with the Fe ioniza- 
tion balance. We computed the weighted mean value of the four 
estimates (ionization balance and Ca lines) as the final estimate 
of the surface gravity logg = 4.66 ± 0.05 dex. 

The metallicity estimate is the mean value of the abundances 
of metallic species with at least 10 spectral lines (Si, Ti, Cr, Fe, 
and Ni); the relative values with respect to the Sun are shown in 
Fig. [5] We note that errors in the element abundance take into 
account the uncertainty in both temperature and gravity, as well 
as the internal scatter of abundances determined from lines of 
the same element. 

The most robust estimates of the effective temperature, sur- 
face gravity, and metallicity are = 5300 ± 60 K, logg = 
4.66 ± 0.09 dex, and [M/H] = +0.39 ± 0.06 dex. For these 
three parameters, systematic uncertainties of 50 K, 0.08 dex, and 
0.05 dex were added. In particular, 40 K was sub tracted from 
r e ff, to be consistent with the temperature scale of iBruntt et al.l 
d2010l) . Finally, the projected rotational velocity was measured 
by fitting synthetic profiles to isolated lines to be vsin/ 
2.0 + 0.8 km s" 1 . This value is higher than the velocity estimated 
from lValenti & Fischerl (120051) . v sin i — 0.9 + 0.5 km s 1 , but it 
still fits inside the error bars. Since both values were extracted 
from spectra of similar resolution (70 000 for Valenti & Fischer 
and 65 000 for NARVAL), we consider the weighted mean value 
v sin i = 1.2 ± 0.4 km s _1 . 

The comparison of the HD 46375 parameters with previous 
wo rks shows reason able good agreement. From spectral analy- 
sis, lLaws et al l d2003l) determined T eff = 5241 ± 44 K, logg = 
4.41 ± 0.09 dex, and [Fe/H] = +0.30 ± 0.03 dex. Moreover, 
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Fig. 4. Comparison of the observed spectrum (grey line) with a 
synthetic profile (blue dashed-dotted line). 
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Fig. 5. Relative abundances of 12 elements measured in 
HD 46375. Open circles and solid squares are mean values for 
neutral and singly ionized lines, respectively. The yellow hor- 
izontal bar marks the mean metallicity with 1-cr uncertainty 
range. 



ISantos et al.ld2004l) found = 5268 ± 65 K, logg = 4.41 ±0.16 
dex, and [Fe/H] = +0.20±0.06 dex. There is good agreement for 
the effective temperature, while our surface gravity and metallic- 
ity are slightly higher than the previous estimates. 

2.4. Stellar luminosity 

From the Hipparcos database, HD 46375 has an a pparent magni- 
tude TO hip = 8.0512 ± 0.0017. IVan Leeuwenl d2007l) reviewed the 
Hipparcos parallaxes and found that n = 28.72±0.89 mas, which 
corresponds to an absolute magnitude of Mhi p = 5.34 ± 0.07 as- 
suming that interstellar reddening is negligible. 

We used b olometric correction s, calculated for the Hipparcos 
spectral band dCavrel et al.l [l997). with the new values of T e g, 
\ogg, and [Fe/H]. The bolometric correction is equal to BC = 
-0.30 ± 0.03. Therefore, the bolometric magnitude is Mboi = 
5.04 ± 0.10, which corresponds to a luminosity L = (0.68 ± 
O.O6)L . 

3. CoRoT observations 

3. 1 . Cleaning of the time series 

As for the previous CoRoT observation runs, the time series at 
level "N2" are still contaminated b y the South Atlantic A nomaly 
(SAA) periodic interruptions (e.g. lAuvergne et al]|2009l) . To di- 
minish the SAA effect, we removed the suspect data and re- 
interpolated all of the gaps by replacing them with an iterative 
estimate of the low f requency trend of the time series, applied by 
iGaulme et alJ d2008l) to Jovian seismic data. 

Since these methods are not able to completely cancel the 
signature of the satellite orbit in the power density spectrum of 
the time series, we eliminate iteratively its harmonics as well 
as the sub-harm onics associated with t he day frequency (11.57 
//Hz) as done in Moss er et a D (|2009b). We cleaned the power 
spectrum across the whole frequency range for every multiple 
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Fig. 6. CoRoT time series of HD 46375 after the re-processing 
described in Sect. 13.11 The time series is plotted in grey at full 
resolution, and in black after a smoothing with box-car of 100- 
min long. The standard deviation around the low frequency trend 
is 176 ppm. The epoch zero corresponds to 3203.62037 days in 
the CoRoT calendar, that is to say after January 1, 2000 at 12:00 
GMT. 



of the satellite frequency (161.72 //Hz) and cleaned 5 daily har- 
monics around each harmonic. 



Density spectrum 

- smoothed density 

- Background fitting 

semi lorentzian components 




10 10" 
v (/iHz) 



Fig. 7. Density spectrum, in log-scale axis. The black curve cor- 
responds to the averaged spectrum obtained with a box-car with 
a moving weighted average, of variable width. There is no clear 
excess power at high frequency that could be directly related to 
pressure modes. 



3.2. Analysis of the stellar activity 

The CoRoT light curve (Fig. |6]l does not exhibit a high activ- 
ity signature, which is also visible in the log-log representation 
of the power density spectrum (Fig. |7). It differs from previ- 
ously observed F-type main-sequence stars (Bar ban et al.l l2009: 



iGarcfa et al.l l2009). Compared to them, there is no clear signa- 
ture at very low fre quencies of t h e rota tion period. A spot anal- 
ysis performed in iMosser et al.l d2009al) is impossible for HD 
46375. The observation time span is not long enough compared 
to the rotation period, and the activity level is very low, as ex- 
pected for this K0V star. It is therefore impossible to identify an 
unambiguous signature of transiting spots. Even after a rebin- 
ning of the light curve at 1 point each CoRoT orbit, the SNR 
remains too low to derive a reliable rotation period. 

Furthermore, the fitting of the power spectrum with a sum 
of 3 semi-Lorentzian profiles, as suggested by Harvey (1985), 
does not converge reliably. Two components are enough to pro- 
vide an acceptable fit. We fitted the power density spectrum with 
a maximum likelihood estimator (MLE), by taking into account 
the nature of the statistics of the power spectrum, which is x 2 
with 2 degrees of freedom. The result is presented in Table 1. 
The amplitude of the power's main component is estimated to 
be 1 . 12*igZ x 1 4 ppm 2 yifHz~ 1 . We note that the upper error bar 
is larger than the value of the amplitude; a longer run and con- 
sequently a higher frequency resolution would allow us to place 
tighter constraints on this parameter. The only reliable estimates 
are the slope pi of each semi-Lorentzian in the log-log represen- 
tation. 



Table 2. Fitted background parameters, where A,, pi, and B 
are defined by the function f(v) = B + A,-/[l + (QiV) Pl ]. 



Parameter 


Value 


Error 


A, 


1.12 10 4 


(+1.67/-0.67) 10 4 


(ppm 2 yiiHz~') 


1.16 10 2 


(+1.97/-0.73) 10 2 


e, 


1.31 10 b 


(+0.51/-0.36) 10" 


(s) 


4.02 10 5 


(+6.3S/-2.46) 10 5 


Pi 


3.94 


±0.30 




1.18 


±0.05 


B (ppmVHz -1 ) 


0.91 


±0.011 



4. p-mode oscillations parameters 

4.1. Estimate of the large separation 

The Fourier spectrum of HD 46375 does not exhibit any seis- 
mic signature visible by eye (Fig. [7]). To determine whether 
th e star displays solar-like oscilla tions, we applied the method 
of Mosser & Appourchaux ( 2009), based on the envelope au- 
tocorrelation function (EACF). This method searches for an 
oscillating signal in the autocorrelation of the time series, 
by analyzing the windowed Fourier spectrum proposed by 
iRoxburgh & Vorontsovl (I2006I) . It has shown to be efficient 
in low SNR cases, such as the CoRoT target HD 175726 
dMosser et al.l l2009bh . The reliability of the result is given by 
an HO test: when the EACF infers a signal above the threshold 
level, the null hypothesis can be rejected. 

The method allows us, in a blind analysis without any a priori 
information, to derive the mean large separation of a solar-like 
spectrum. We were then able to derive an unambiguous detec- 
tion around 155 pHz. The autocorrelation factor, around 1 1.5, is 
much above the threshold level of 8 calculated for a blind detec- 
tion. 
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Table 3. Classical and seismic parameters of HD 46375 derived 
from the observations and the model that minimizes the^ 2 func- 
tion defined in Sect. 5.3. 



Fig. 8. Envelope of the autocorrelation function as function of 
the frequency and the time. The maximum value region, in the 
frequency range [2600, 3400] yuHz, corresponds to a correlation 
time of about 3.6 h, that is to say a large separation of 153 
//Hz. The EACF value is larger than 8 in the frequency range 
[2750, 3300] //Hz. 



2.6 2.8 3.0 3.2 3.4 3.6 

v (mHz) 

Fig. 9. Variation in the large separation with frequency, derived 
from the autocorrelation of the time series. The power spec- 
trum was first windowed with a 0.75 mHz-half width cosine 
filter, as indicated in the lower-right corner. The sampling of 
the curve corresponds crudely to the mean large separation. The 
size of the dots is representative of the amplitude of the corre- 
lated signal. One-cr error bars (dashed lines) are derived from 
Mosser & Ap pourchauxl (120091) . 



4.2. Variation of the large separation with frequency 

Measuring eigenfrequencies in the Fourier spectrum of 
HD 46375 remains impossible, even by folding and collapsing 
the po wer spectrum in a collapsogram, to display a comb-like 
signal (Korzennik 2008). However, the autocorrelation method 
makes it possible to investigate the variation Av(v) in the large 
separation with frequency. The performance of the method de- 
pends on the autocorrelation factor. In the case of HD 46375, 
we can obtain 2 independent measures of the large separation, 
as indicated by Fig. [8] from which we derive the large separa- 
tion increase with frequency. The analysis of the spectrum also 
provides error bars, as indicated in Fig. [9] 



4.3. p-mode amplitude 

The maximum bol ometric amp l itude of radial modes was de- 
rived according to Mic hel et al.1 (120081) . The amplitude appears 



Parameters 


Observations 


Model 


Observational constraints 


T eS (K) 


5300 ± 60 


5326 


log(L/L e ) 


-0.17 ±0.04 


-0.22 


l°gg 


4.66 ± 0.09 


4.51 


(Z/X\ 


0.060 ± 0.008 


0.0626 


<Av> 


153.0 ± 0.7 


153.4 


Free parameters 


Mass (M Q ) 




0.97 ± 0.05 


Age (Gyr) 




2.6 ±0.8 


Y 




0.31 ±0.05 


x 1 




2.51 



to be quite low, about 1.25 + 0.25 ppm. The large uncertainty 
originates first from the difficulty to estimate the background, 
and second from the very low height-to-background ratio %, of 
about 2.5 %. This ratio measures the oscillation power relative to 
the power of the background, both intergrated over a 2-Av wide 
frequency range centered around the maximum oscillation fre- 
quency. The actual value of 2.5 % is very low, indeed much lower 
than the minimum value for which solar-like oscillations were 
derived in pr evious CoRoT obser vations (0.11 for HD 175726 
according to Moss er et al.l 1200 9b). This shows the difficulty in 
measuring the large separation of HD 46375. 

We can compare the ob served am plitudes to the scaling law 
(L/M) 0J proposed bv lSamadi et all (120071) for the Doppler sig- 
nal. When translate d into a photometric amp litude thanks to the 
adiabatic relation (Kield sen & Beddirial 19951) . this gives a scal- 
ing low (L/M) 0J I -yTeff according to which the bolometric am- 
plitude of radial modes are estimated at the level of A\, \(£ = 
0) = 1 .7 ppm. As in most of the main-sequence stars observed by 
CoRoT, the observed amplitudes are significantly below the pre- 
dictions. Mode lifetimes cannot be measured. By comparing the 
frequency amplitudes with the autocorrelation signal, one can in- 
fer that mode lifetimes are longer at low frequency, as observed 
in other stars. 



5. Modelling 

We now construct a detailed asteroseismic model of HD 46375. 
The aim is to identify which stellar models are consistent with 
both the classical fundamental parameters (T e ff , log g, and metal- 
licity) from Sect. I2.3l and the seismic constraints from Sect. [4] 

5.1. Seismic constraints 

The autocorrelation method enabled us to obtain an estimate of 
the variation in the large separation with frequency. However, the 
frequencies at which we estimated Av depend on the width of the 
filter that we used, and the large separations cannot be directly 
compared to those of a model. We consequently used a mean 
value of the large separation (Av) over the frequency domain 
where the EACF produces a signal above the threshold level. We 
also performed a linear regression upon Av, taking into account 
the corresponding uncertainties given in Fig. [9] To perform a 
rapid but r obust modeling, we reduced this variation to a linear 
trend, as in Mosser et al. (2008). We therefore consider (Tabled 



Av = Avo + y(v - vo) 



(3) 
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The frequency vo = 3 mHz corresponds to the frequency of the 
maximum amplitude. 

5.2. Properties of the models 

Models were computed with the evolution code CESAM2k 
dMorell 1 19971) . and the mode frequencies were derived from 
the models using th e Liege OScillation Code (LOSC, see 
IScuflaire et all 120081) . We used the O PAL 2001 equation o f 
state and opacity tables as described in Lebretonet al. (2008). 
Nuclear reactions rates were computed using the NACRE com- 
pilation (Angulo et al. 1999). The models were calcu lated with 
the element abundances of iGrevesse & Noelsl 0993!). We used 
Eddington's grey temperature-optical depth law to describe the 
atmosphere. The connection with the envelope was achieved 
at an optical depth of 10. Convection was described using the 
Canuto, Goldman, and Mazzitelli (CGM, ICanuto et all 1 19961) 
formalism, which involves a free parameter assumed to be some 
fraction cccgm of the local pressure scale height H p . In our study, 
this parameter was fixed to the solar value (q'cgm = 0.64). 

Overshooting is described as an extension of the convective 
core over a distance d ov = a m min(r cme ,Hp), where r core is the 
radius of the convective core. Microscopic diffusion is not taken 
into account. 

In contrast to heavier elements, it is difficult to determine 
spectroscopically the helium content F of the stellar surface. 
Hence, Y is assumed to be a free parameter of the models. 
However, several studies have shown that a link exists between 
the helium enrichment of the galaxy and the enrichment of heav- 
ier metals. This leads to the determination of a helium-to-metal 
enrichment ratio, written AF/AZ such that 



AY 

Y = Y P + — Z. 
AZ 



(4) 



For K-type stars ar ound and above solar metallicity, 
ICasagrande et all d2007l) found AF/AZ ~ 2.1 ± 0.9. We can thus 
derive an estimate of F. Inserting Eq.|4]into X + Y + Z = 1, one 
indeed obtains 



Y= 1- 



1-Fp 



1 + (Z/X) (1 + AF/AZ) 



[1 + (Z/X)] 



(5) 



Taking F P = 0.248 ± 0.003 (see iPeimbert et aTll2007l) . and us- 
ing our spectroscopic determination of the metallicity (Z/X) = 
0.063 ± 0.006, we find that F = 0.33 + 0.05. 

5.3. Optimization 

We now search for models that reproduce the observational con- 
straints (the effective temperature T s g, the luminosity L, the sur- 
face gravity \ogg, the metallicity (Z/x) s , and the mean large sep- 
aration (Ay)) using as free parameters the stellar mass, age, and 
initial content in helium Yq. The values of the observational con- 
straints are summarized in Table 3. 

To estimate the compatibility of the models with the obser- 
vations, we use 



X 



Zl 



A° bs - A 



mod \2 



(6) 



where A° hs is the i th observational constraint (among N), cr° hs its 
l-cr error bar, and /t™ od the value of the corresponding parameter 
in the considered model. 



To search for an optimal model, we use the Levenberg- 
Marquardt algorithm, as described by Mi ali o & MontalMrJ 
d2005l) . The interest of this method is that it enables us to min- 
imize a function with a small number of iterations, even with 
initial parameters rather far away from the final result. However, 
when using this optimization method, we must pay attention to 
the risk of converging to a local, not global, minimum. Hence, 
we combined this approach with the computation of a grid of 
models. 

The parameters of the optimal model we found are given in 
Table 3. All the classical parameters are fitted within l-cr er- 
ror bars, except the stellar luminosity for which the obtained 
value differs from the observed one by 1 .4 cr. The mean value 
of the large separation is reproduced well by our best-fit model. 
We however note that if we perform a linear fit to the varia- 
tions in the large separation for this model, we obtain a slope 
of 7mod = —1.1 10~ 4 , quite different from the one we derive 
with the autocorrelation of the time series (0.010 + 0.003). To 
investigate this disagreement, individual frequencies would be 
required. This would require more complex modelling based on 
observations in higher SNR, as should be provided with the re- 
observation of this exoplanet-hosting star. 

6. Discussion and conclusion 

This work is particularly interesting because of its overlap 
between exoplanetary science an d asteroseismology, as also 
perfo rmed by the Kepler mission dChristensen-Dalsgaard et all 
2010) and the proposed Plato mission. Even for very low signal- 
to-noise ratio data, as presented here, we have demonstrated that 
we can constrain the stellar fundamental parameters and hence 
also the planetary parameters. Moreover, the ground-based sup- 
port data provided additional constraints on the magnetic field 
as do high resolution spectra, which have been used to refine the 
atmospheric parameters. 

From the Zeeman Doppler imaging observations, performed 
with the NARVAL spectropolarimeter at the Pic du Midi 
Observatory, we reconstructed the stellar magnetic structure. 
We found that the magnetic field of HD 46375 is signifi- 
cantly different from other planet-bearing stars, such as t Boo 
and HD 189733: for HD 46375, the magnetic field is mainly 
poloidal, while it is predominantly toroidal for the two other 
stars. This difference is related to the stellar rotation. Using the 
rotation-activity relationship, the stellar rotation period is 42^ 
day. Moreover, since the magnetic field measurements do not 
display any periodicity, we have confirmed that the period is 
longer than the run-duration (35 days). With this rotation period, 
HD 46375 is among stars featuring very simple, mostly poloidal, 
and axisymmetric large-scale magnetic topologies. 

The comparison of the fundamental parameters before and 
after this work illustrates the potential of asteroseismology. 
The stellar mass has increased by about 7 % and the stel- 
lar radius by the same proportion, which corresponds to a 
much denser star with a correspondingly higher large separa- 
tion. Desp ite the low signal, the runnin g filtered autocorrelation 
method (Mosser & Aprjourchaux 2009) allows us to measure the 
large separation of the p-mode frequencies. The mean value of 
the large separation is 153.0 ± 0.7 //Hz, which is the highest 
value of the large separation ever observed in photometry. Only 
three targets with higher Av have been observed in spectroscopy 
(t Ceti, a Cen B, and 70 Oph A). The Fourier analysis of the 
CoRoT light curve exhibits a tiny excess of power, with am- 
plitude of about 1.25 ppm, i.e., a height-to-background ratio of 
about 2.5%. 
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We combine seismic parameters with the atmospheric pa- 
rameters obtained from the high resolution NARVAL spectra. 
This allows us to refine the fundamental parameters of the star 
i.e., its mass, effective temperature, radius, and age. The best-fit 
model has a mass of 0.97 ± 0.04 Mq with an age of 2.6 ± 0.8 
Gyr. Despite the relatively large error bar in the age, all the other 
inputs, such as the high metallicity, the lower temperature, and 
luminosity with respect to the Sun reinforce the hypothesis of 
a young metal-rich star. Thus, it is not exactly a solar twin, 
but a young Sun with a composition enhanced in metallic ele- 
ments. As expected from the a steroseismic diagram proposed by 
Christensen-Dalsgaard (1988), this younger analog of the Sun 
has a higher value of the large separation. However, since we 
have not been able to identify the individual eigenfrequencies, 
we cannot constrain the age more tightly. 

The revision of the stellar mass permits us to re-estimate 
the projected mass of the known exoplanet to b e M D sin ; = 
0.234 + 0.008 Mj up , instead of 0.226 ±0.019 M jup in lButler et all 
(2006). Although we claim an improvement in the accuracy by a 
factor two, we cannot constrain the true planetary mass. If we as- 
sume the planet orbital plane is in the stellar equatorial plane, the 
planetary mass is M p = 0.29!°;^ Mj U p. However, we do not know 
whether the orbital plane is in the stellar equatorial plane as in 
the Solar System: seve ral exo-planetary systems do show diff er- 
ent configurations fe.g. lWinn et alJ2009t iSibrard et alJ[2008h . 

The CoRoT short run on HD 46375 was part of an additional 
program: it has demonstrated that we are able to carry out a com- 
plete analysis of a stellar system when combining photometry 
and ground-based spectrometric support. However, a 34-day run 
appears too short to provide high-quality data for such a rela- 
tively faint star. From a 5 -month long run with CoRoT, it would 
be possible to derive individual eigenfrequencies and to monitor 
in more detail the rotation period. Moreover, a longer run should 
be able to confirm the detection of a planetary contribution to the 
light curve, to be reported in our forthcoming Paper EL 
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